Propagation of Avalanches in Mni 2 -acetate: Magnetic Deflagration 
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Local time-resolved measurements of fast reversal of the magnetization of single crystals of Mni2- 
acetate indicate that the magnetization avalanche spreads as a narrow interface that propagates 
through the crystal at a constant velocity that is roughly two orders of magnitude smaller than the 
speed of sound. We argue that this phenomenon is closely analogous to the propagation of a flame 
front (deflagration) through a flammable chemical substance. 



Mni2-acetate (hereafter Mni2-ac) is a prototypical 
molecular magnet composed of magnetic molecules, 

[Mni20i2(CH 3 COO)i 6 (H20)4]-2CH 3 COOH4H20, 
with cores consisting of twelve Mn atoms strongly cou- 
pled by exchange to form superparamagnetic clusters of 
spin S — 10 at low temperatures-. Arranged in a cen- 
tered tetragonal lattice, the spin of the Mni2 clusters is 
subject to strong magnetic anisotropy along the symme- 
try axis (the c-axis of the crystal). Below the blocking 
temperature of ~ 3.5 K, the crystal exhibits remarkable 
staircase magnetic hysteresis due to resonant quantum 
spin tunneling between energy levels on opposite sides 
of the anisotropy barrier corresponding to different spin 
projections, as illustrated in Fig. l(a)£ This and other 
interesting properties of Mni2-ac have been intensively 
studied in the last decade (see Refs. ytllalafor reviews). 

It has been known for some time! that Mni2-ac crys- 
tals exhibit an abrupt reversal of their magnetic moment 
under certain conditions. This phenomenon, also ob- 
served in other molecular magnets, has been attributed 
to a thermal runaway (avalanche) in which the initial 
relaxation of the magnetization toward the direction of 
the field results in the release of heat that further accel- 
erates the magnetic relaxation. Direct measurements of 
the heat emitted by Mni2-ac crystals^ as well as mea- 
surements of the magnetization reversal in pulsed mag- 
netic fields^ have confirmed the thermal nature of the 
avalanches. More recently, the electromagnetic signal as- 
sociated with avalanches was detectedi2*ii and it was ar- 
gued that if the radiation is of thermal origin it would 
indicate a significant increase in the temperature of the 
crystal. This has not been confirmed by direct bulk mea- 
surements of the temperature using a thermometer. Ev- 
idence has been obtained^ that the avalanche may not 
be a uniform process throughout the sample. No clear 
understanding of the avalanche process has emerged to 
date. 

In this Letter we report local time-resolved measure- 
ments of fast magnetization reversal (avalanches) in mm- 
size single crystals of Mni2-ac. We show that a mag- 
netic avalanche takes the form of a thin interface be- 
tween regions of opposite magnetization which prop- 
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FIG. 1: (a) The potential energy function in a longitudinal 
magnetic field; (b) Schematic diagram of a sample mounted 
on an array of Hall sensors used to detect B x . The mag- 
netization is relatively uniform for quasi-equilibrium condi- 
tions; pronounced non-uniformity in magnetization during an 
avalanche generates large values of B x . 



agates throughout the crystal with a constant field- 
dependent speed ranging from 1 to 15 m/s. We demon- 
strate that this phenomenon is closely analogous to the 
propagation of a flame front (deflagration) through a 
flammable chemical substance. 

Microscopic arrays of Hall bars were used to measure 
the magnetization of three single crystals of Mni2-ac with 
dimensions: sample 1, 0.29 x 0.29 x 0.64 mm 3 ; sample 2, 
0.28 x 0.28 x 1.44 mm 3 ; sample 3, 0.24 x 0.24 x 1.02 mm 3 . 
Eleven Hall bars of dimensions 10 x 10 /im 2 with 30 /j,m 
intervals were used for sample 1, and 30 x 30 /zm 2 with 
130 /im intervals for samples 2 and 3. 

Using an excitation current of 2 /iA, the Hall bar sig- 
nal was amplified by a factor of 1200, and detected and 
recorded by several digital scopes and a data acquisi- 
tion card. In order to ensure proper synchronization, one 
channel of each scope was anchored to the same signal. 
The Hall sensor and amplifier introduced combined de- 
lays of up to 3 [is. A magnetic field was applied in the z- 
direction along the crystal easy axis (see Fig. 1), lowering 
(raising) the energy of the states corresponding to spin 
projections along (opposite to) the field direction. The 
Hall sensors were aligned to detect the magnetic induc- 
tion of the sample in the ^-direction. B x is proportional 
to the spatial derivative of M z in the region near the sen- 
sor. For a uniform magnetization in the z-direction, B x 
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FIG. 2: Hysteresis loops for a single crystal of Mni2-ac, and 
an avalanche; data shown for sample 1. B x was recorded by 
a Hall sensor situated about 80 /im above the middle of the 
sample. Steps due to quantum tunneling of the magnetization 
are observed, followed by a sharp spike in B x associated with 
an avalanche, shown in the inset. 




FIG. 3: Signals recorded by seven equally spaced Hall sen- 
sors situated near the center of sample 1 during an avalanche 
triggered at 4 T. The left (right) trace corresponds to the top 
(bottom) sensor (see Fig. 1). Sensor positions are measured 
relative to the top of the sample, with the center at ~ 320 
fim. The inset shows the sensor position versus the time at 
which the sensor recorded peak amplitude. A straight line fit 
yields a constant velocity of propagation of 12 m/s. 



derives from the gradient at the sample ends, which is 
proportional to M z itself. During an avalanche, there is 
a large contribution to B x from the local region corre- 
sponding to the avalanche front, where dM z /dz is large. 

The samples were immersed in liquid 3 He; most of the 
data were obtained at the base temperature of 250 mK. 
The few points measured at 400 and 650 mK were found 
to lie on the same curve within the scatter of the data, 
indicating that the temperature dependence is weak. A 
longitudinal magnetic field (parallel to the easy axis) was 
swept back and forth through the hysteresis loop to ±6 
T until an avalanche was triggered. As reported in an 
earlier paper—, avalanches occur in a stochastic way at 
0.25 K both at resonant magnetic fields (where energy 
levels on opposite side of the barrier match, see Fig. 1) 
and away from resonance. Avalanches were also found 
for sample 2 for zero field-cooled conditions, where the 
sample starts from zero magnetization (instead of full 
saturation) . 

Figure 2 shows an avalanche for sample 1. Steps due to 
quantum tunneling of the magnetization were observed, 
with a magnetization that was almost uniform through- 
out the sample, until an avalanche occured, as shown in 
the inset. During the avalanche the Hall bar recorded a 
large peak in B x , signaling the abrupt onset of a highly 
non-uniform magnetization. 

For a field sweep rate of 10 mT/s and temperature 
0.25 K, Fig. 3 shows an avalanche triggered at 4 T 
and recorded for sample 1 by seven of the eleven sen- 
sors placed in sequential positions near the center of the 
sample. The avalanche was triggered above the top-most 
sensor and traveled downward (see Fig. 1). B x displays 
the largest peak at the center due to the finite size of the 
sample. The inset shows the sensor position as a func- 
tion of the time at which the sensor registered the peak 




HoH(T) 

FIG. 4: Velocity of propagation of avalanches versus longi- 
tudinal magnetic field at which the avalanches occured. ZFC 
denotes data obtained for sample 2 cooled in zero field, thus 
starting from zero magnetization and ending in full magneti- 
zation, AM/2M aa t = 0.5; AM/2M sat = 1 for sample 3 and 
the remaining sample 2 data; for sample 4, AM /2M sa t varied 
around 0.7. From top to bottom, the curves are fits to the 
data for AM/2M sat = 1,0.7 and 0.5 (see text). 



amplitude. The slope of the straight line drawn through 
these points yields a constant velocity of 12 m/s for this 
avalanche. 

Figure 4 summarizes the data obtained for the velocity 
of propagation of avalanches recorded in different longitu- 
dinal magnetic fields for all three samples. For avalanches 
starting from full magnetization, the data for samples 2 
and 3 lie on approximately the same curve. The velocity 
decreases with decreasing longitudinal magnetic field and 
goes to zero at about 0.6 T, below which no avalanches 
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through which the heat diffuses during the time of the 
"chemical reaction" r. In our case 
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FIG. 5: Velocity of propagation of avalanches versus energy 
released per molecule. The line is drawn to guide the eye. 



can occur. Smaller velocities are obtained for avalanches 
in sample 2 when starting from the zero-field-cooled con- 
dition. Avalanches for sample 1 were obtained only at 
relatively high magnetic fields in the vicinity of 4 T; the 
velocities for this sample range in value and do not appear 
to be consistent with data for the other two samples. 

Interestingly, as shown in Fig. 5, an approximate col- 
lapse is obtained for all the data when plotted as a func- 
tion of g^BHS(AM/M sat ), the energy per molecule re- 
leased during an avalanche. Thus, avalanches require the 
release of a threshold energy, above which they propagate 
with a speed that appears to be a linear function of the 
energy for the range investigated in these experiments. 

Our observations cannot be attributed to magnetiza- 
tion reversal associated with domain wall motion, since 
there is no long-range order in our system. Some insight 
can be obtained by noting that, from a thermodynamic 
point of view, a crystal of Mni2 molecules placed in a 
magnetic field opposite to the magnetic moment is equiv- 
alent to a metastable (flammable) chemical substance. 
In our case, the role of the chemical energy stored in a 
molecule is played by the difference in the Zeeman en- 
ergy, AE = 2g/iBHS, for states of the Mni 2 -ac molecule 
that correspond to S parallel and antiparallel to H; here 
g = 1.94 is the gyromagnetic factor and is the Bohr 
magneton. For Mni2-ac in a field of a few Tesla, AE 
is below 0.01 eV, as is the energy barrier, U(H), be- 
tween spin-up and spin-down states due to the magnetic 
anisotropy. Thus, for the avalanches in Mni2-ac, AE and 
U are two orders of magnitude smaller than typical en- 
ergies of chemical reactions. However, our temperature 
range is also more than two orders of magnitude below 
room temperature, making the analogy rather close. 

A well-known mechanism for the release of energy by 
a metastable chemical substance is combustion or slow 
burning, technically referred to as deflagrationiA 5 It oc- 
curs as a flame of finite width, 6, propagates at a constant 
speed, v, small compared to the speed of sound. The 
parameter S is determined by the distance, S ~ \[kt 



t(H) = t exp 



k B T f 



(1) 



where tq ~ 10~ 7 s is the attempt time^ and Tf is the 
temperature of the flame. The dynamics of the flame are 
governed by the thermal diffusivity, k, which obeys: 



(2) 



For k independent of T, substituting T — T(x — vt) at 
x > vt, one obtains T — Tf exp [—v(x — vt)/ k] in front of 
the interface, which yields vS — k. An interface thickness 
that is at most the distance between sensors, 5 ~ 30/xm, 
and the experimentally measured velocities of the or- 
der of 1 — 15 m/s, yield an upper bound on k in the 
range 10 -5 m 2 /s to 10~ 4 m 2 /s, consistent with heat pulse 
experiments^ 

Combining vS — k with 6 ~ \/kt , one obtains: 



1/2 



exp 



U(H) 



2k B T. 



(3) 



The strongest dependence of v on H derives from the 
exponential, which contains the known dependence^ of 
the energy barrier, U(H), on the magnetic field. Assum- 
ing that the temperature of the flame is proportional to 
the released magnetic energy density, Tf = CHAM, it is 
possible to fit all the data with a single value of the pro- 
portionality constant C, as shown in Fig. 4. The flame 
temperature obtained from these fits ranges from 8.5 K 
for an avalanche at 1 T to 26 K for an avalanche triggered 
at 4 T. We note that the prefactor obtained from the fit is 
consistent with the values of k and To discussed earlier^ 
Here we have used the simplest model of deflagration, 
a widely studied phenomenon that is known to be quite 
complex^ This crude model captures the overall behav- 
ior, and yields parameters that are quite reasonable in 
size. A more complete theory is needed to account for 
the apparent data collapse of Fig, 5. 

Recent bolometer measurements of the radiation gen- 
erated by a magnetic avalanche 10 ' — gave puzzling re- 
sults that can be understood within our model. One 
enigma was that the sample temperature measured by a 
thermometer directly following an avalanche was lower 
(< 6 K) than the temperature registered by the bolome- 
ter if one assumed thermal radiation. A second puzzle 
was that the reversal of the magnetization during the 
avalanche occured on a much shorter time scale than 
the cooling of the sample following the avalanche>ii We 
suggest that the radiation observed during avalanches is 
generated by the narrow hot interface (flame) that prop- 
agates through the crystal. The temperature of the bulk 
of the crystal (including the "ash" left behind the inter- 
face) is always significantly lower than the temperature 
of the interface itself. The time of the magnetization re- 
versal is determined by the time t = l/v needed for the 
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interface to sweep the sample of length I. In our case 
t ~ 0.1ms, while the time needed for the "ash" to reach 
equilibrium with the thermal bath can be much longer. 

The strongest evidence that our observations are due 
to deflagration is the presence of a well defined propa- 
gating front requiring a threshold energy traveling at a 
subsonic velocity. The deflagration mechanism provides 
the condition needed for the avalanche to occur. This 
condition is the same as the condition needed to sustain 
the propagation of a flame through a chemical substance. 
It is well-known that deflagration of a flammable gas will 
not occur in a pipe of diameter, d, small compared to 
the width of the flame, S. In our case S must be small 
compared to the diameter of the crystal. If this condition 
is not satisfied, the heat generated by the magnetization 
reversal diffuses mostly through the walls of the sam- 
ple and cannot sustain the propagation of the interface. 
This explains why avalanches only occur in larger crys- 
tals with sufficiently large magnetization opposite to the 
direction of the field. The latter condition coincides with 
the condition of "flammability"— needed to provide suf- 
ficient heating (that is, the large Tf) required for 5 < d. 
It is interesting to note in this connection that the few 
avalanches that were recorded around 1 T did not result 
in full reversal of the magnetization. At these low fields 
the conditions for ignition are marginally satisfied, and 
the "flame" is extinguished before the process of mag- 
netization reversal has been completed. In addition to 
available magnetic energy (flammability), the conditions 
for ignition may also depend on the shape and quality of 
the crystal, which may account for differences observed 
for different samples. 

Slow burning at a subsonic speed (deflagration) is gov- 



erned by the linear process of thermal conductivity. In 
addition to deflagration, unstable chemical substances 
also exhibit detonation, which can be caused by instabil- 
ity of the flame or by direct initiation other than through 
deflagration*^ The initial stage of the detonation corre- 
sponds to a non-linear supersonic shock wave>^*i^ The- 
ory and experimental studies of advanced stages of det- 
onation are lacking. Based on the close analogy between 
unstable chemical substances and molecular magnets, the 
latter may well exhibit "magnetic detonation" under the 
right conditions. 

In conclusion, we have demonstrated that avalanches 
in the magnetization reversal of sufficiently large crystals 
of magnetic molecules are very similar to flame propa- 
gation (deflagration) through a metastable chemical sub- 
stance. The analogy between the two systems derives 
from the magnetic bi-stability of molecular nanomag- 
nets. Our observation of "magnetic deflagration" offers a 
potentially important new way to investigate the phe- 
nomenon of flame propagation (and, possibly detona- 
tion). In contrast to deflagration in flammable chemi- 
cal substances, the analogous process of "magnetic de- 
flagration" in molecular nanomagnets is non-destructive, 
reversible, and much easier to control. 
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The prefactor of the exponential in Eq. [3] the second pa- 
rameter required for the fits, differed by as much as 30% 
for different AM, reflecting the fact that the prefactor is 
not well known within the theory. On the other hand, the 



exponential dependence is well captured by the model. 
I. Glassman, Combustion (Academic Press, f996). 



